Evidence for Xcitt" at Belle 
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We present a measurement of the branching fraction for the Cabibbo- and color-suppressed B" 
Xci7r° decay based on a data sample of 657 x lO'' BB events collected at the T(45') resonance with the 
Belle detector at the KEKB asymmetric-energy e'^e~ collider. We observe a signal of 40 ± 9 events 
with a significance of 4.7o" including systematic uncertainties. The measured branching fraction is 
B{B° Xci7r°) = (1.12 ± 0.25(stat.) ± 0.12(syst.)) x 10"^ 

PACS numbers: 13.25.Hw, 14.40Gx, 14.40.Nd 



The decay B'^ Xcitt^ is a 6 — > ccd transition that 
proceeds at leading order through the color-suppressed 
tree diagram as shown in Fig. 1. If the tree dia- 
gram dominates, then the time-dependent CP-violating 
asymmetries in this decay mode are predicted to be the 
same as those measured in 5 ^ ccs decays, such as 
deviation of the CP-violating asym- 
metries in P° XciT^^ from these expectations could 
indicate non-negligible contributions from a penguin am- 
plitude or amplitudes from new physics. For a similar 
B decay mode, P° J/tpir'^, the time-dependent CP- 
violation parameters have been measur ed by the Belle [3 
and BaBar Q collaborations. Comparison of the prop- 
erties of P° J /ipTT^ and B^ Xcitt^ decays will also 
provide an opportunity to probe new physics that pre- 
dicts different couplings to left-handed and right-handed 
particles Q. 

The P° Xci7i'° decay has not been observed so far. 
Confirming its existence is a very important step toward 
detailed studies of the b — s- ccd transition. The factoriza- 
tion approach [5| and isospin symmetry imply that the 
branching fraction of the P° Xciti"" decay mode should 
be one half of that for B^ Xcitt^ Precise measure- 
ment of the branching fractions of these decays can also 
provide information related to the final state interactions 
in B decays. 

In this paper, we report the first evidence of P° 
Xci7r° using a data sample containing (657±9) x 10^ BB 
events collected at the T(45) resonance with the Belle 
detector f?\ at the KEKB asymmetric-energy e+e~ col- 
lider Q. The Belle detector is a large solid-angle mag- 
netic spectrometer located at the KEKB e+e" storage 
rings, which collide 8.0 GeV electrons with 3.5 GeV 
positrons producing a center-of-mass (CM) energy of 
10.58 GeV, the mass of the T(4S') resonance. 

The Belle detector consists of a silicon vertex detector 
(SVD), which is surrounded by a 50-layer central drift 
chamber (CDC), an array of aerogel Cherenkov counters 
(ACC), a barrel- like arrangement of time-of-flight (TOF) 
scintillation counters, and an electromagnetic calorime- 
ter (ECL) comprised of CsI(TZ) crystals. These subde- 
tectors are located inside a superconducting solenoid coil 




FIG. 1: Leading tree level diagram for B^ —>■ Xci7i"° decay. 



that provides a 1.5 T magnetic field. An iron flux-return 
yoke (KLM) located outside the coil is instrumented to 
detect Kl mesons and to identify muons. The detector 
is described in detail elsewhere \r\. The data set used in 
this analysis consists of two subsets: the first 152 x 10^ 
B meson pairs were collected with a 2.0 cm radius beam 
pipe and a 3-layer SVD, and the remaining 505 x 10^ B 
meson pairs with a 1.5 cm radius beam pipe, a 4-layer 
SVD and a small-cell inner drift chamber [l3| . 

Events with B meson candidates are first selected 
by applying the following general selection criteria for 
hadronic events: at least three charged tracks are re- 
quired to originate from an event vertex which is con- 
sistent with the interaction point (IP); the reconstructed 
CM energy should satisfy iJ"^*^ > 0.2y^, where y/s is the 
total CM energy; the component of momentum along 
the beam direction (z-direction) must be in the range 
Ip^^'^I < 0.5y/s/c; and the total ECL energy should con- 
sist of at least two energy clusters and satisfy 0.1 ^/s < 
^ECL < 0-8y/s. To suppress continuum background 
dominated by two-jet-like e~^e~ qq annihilation (q = 
u, d, s), we reject events where the ratio of the second to 
zeroth Fox- Wolfram moments [ll| i?2 is greater than 0.5. 
We find no contribution from continuum background af- 
ter applying this cut on R2 . To remove tracks of charged 
particles that are poorly measured or do not originate 
from the interaction region, we require their origin to be 
within 0.5 cm of the IP in the radial direction, and 5 cm 
in ^-direction. 

We reconstruct J/tp from pairs of e^e~ or can- 
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FIG. 2: The AM (Mi+f-^- Mf+e-) distribution for inclusive 
Xci candidates. The enhancement just above the Xci mass 
region is due to the Xc2- 



didates. For muon tracks, identification is based on track 
penetration depth and the hit pattern in the KLM sys- 
tem. Electron tracks are identified by a combination of 
dE/dx from the CDC, E/p [E is the energy deposited 
in the ECL and p is the momentum measured by the 
SVD and the CDC), hght yield in the ACC, shower 
shape in the ECL and position matching between an ECL 
cluster and extrapolated track. In order to recover di- 
electron events in which one or both electrons radiate a 
photon, the four-momenta of all photons that lie within 
0.05 radians of the e"*" or e~ directions are included in 
the invariant mass calculation. The invariant mass win- 
dow used to select J/i/j candidates in the /j,"'"^~(e"'"e^) 
channel is -0.06 (-0.15) GeV/c^ < M^+^- - ruj/^ < 
0.036 GeV/c^, where nij/^ denotes the nominal J/ip 
mass 12]; these intervals are asymmetric in order to 
include part of the radiative tails. Vertex- and mass- 
constrained kinematic fits are performed for selected J/tp 
candidates to improve the momentum resolution. 

Photons are identified as ECL energy clusters that are 
not associated with a charged track and have a minimum 
energy of 0.060 GeV. We reject the photon candidate if 
the ratio of the energy deposited in the array of the cen- 
tral 3x3 ECL cells to that in the array of 5 x 5 cells 
is less than 0.82. Neutral pion candidates are formed 
from the photon pairs that have an invariant mass in the 
range 0.118 GeV/c^ to 0.150 GeV/c^. To reduce com- 
binatorial background, an energy asymmetry threshold 
\{E^i — Ej2)/iE-yi + E^2)\ < 0.8 is used, where Eji and 
E^2 are the energies of each photon of the tt" candidates 
in the CM frame. Finally, a mass constrained fit is ap- 
plied to the 7r° candidates. 

To reconstruct the Xci meson, we combine a J/^p candi- 
date with momentum below 2.0 GeV/c in the CM frame 
with a selected photon. To suppress photons originat- 



ing from TT 77, we veto photons that, when combined 
with another photon in the event, satisfy 0.110 GeV/c^ < 
M^-y < 0.150 GeV/c^. The Xci candidates are selected by 
requiring the mass difference {AM — Mi+i-^ — Mf+i-) 
to lie between 0.3 GeV/c^ and 0.5 GeV/c^. The AM 
distribution is shown in Fig. 2. A mass-constrained fit 
is applied to Xci candidates in order to improve the mo- 
mentum resolution. 

We reconstruct B mesons by combining a Xci can- 
didate with a neutral pion. The energy difference, 
AE = E*g - E*^^^ and the mass difference (AM = 
Mi+i-^ — Mf+i- ) are used to separate signal from back- 
ground, where i?boam ^^'^ run-dependent 
beam energy and reconstructed energy of the B meson 
candidates in the CM frame, respectively. For the se- 
lected 5° candidates, the beam-constrained mass, Mbc = 
V ^bcam ~ P*B ' whcrc p*g IS the reconstructed momentum 
of the B meson candidates in the CM frame, is required 
to be 5.27 GeV/c2 < Mbc < 5.29 GeV/c^. We retain 
candidates with 0.3 GeV/c^ < AM < 0.5 GeV/c^ 
and -0.2 GeV < A£; < 0.2 GeV for the final analysis. 
The selection criteria are determined by optimizing the 
figure of merit, S/y/{S + B), where S (B) is the num- 
ber of signal (background) events in the signal region 
(-0.09 GeV < AE < 0.05 GeV and 0.380 GeV/c^ < 
AM < 0.435 GeV/c^). We assume the signal branching 
fraction is half of that for B^ — > Xcitt^ 0. 

The backgrounds are dominated by the BB events 
with a J/f/; in the final state, where the J/ip is pro- 
duced either directly from B decay or from the Xci ~* 
decay chain. We study these backgrounds using 
a large Monte Carlo (MC) sample [3] corresponding to 
3.86x10^" generic BB decays that includes all known 
B J/ipX processes and those where B decays into 
higher charmonium states (xci, Xc2 or ip') that subse- 
quently produce J/ip in the final state. The dominant 
contribution comes from B° ^ J/ipK°s{^ Tr°Tr°), B± ^ 
J/2PK*^{892), B° J/^PK*°{892), ^ XciK^{-^ 
7r°7r°), and a few other exclusive B — > J/V'(Xci) + X de- 
cay modes. To suppress neutral pions from Kg — > tt^tt*^ 
decays, we veto tt^s that, when combined with another 
7r° in the event, satisfy 0.469 GeV/c^ < M^o^o < 
0.526 GeV/c^. The Kg veto reduces the background 
by 16.2% with a signal loss of 4.2%. We further re- 
duce the B J/4'Kg{-^ 7r°7r°) background by requiring 
cos^hci > —0.88, where 6'hci is the angle between the di- 
rection opposite to the B momentum and the 7 direction 
in the Xci rest frame. After this requirement, we find that 
there is no peaking background in the AE signal region. 
The background from B — > Xci + X decay modes, such 
as B° ^ Xcii^g(^ ttOttO), B± ^ Xci^^*, ^ XciP^ 
and B" -> ip{2S){-^ Xci7)'^°, forms a peak in the AM 
signal region (called peaking background) , while all other 
components are fiat in AM (called combinatorial back- 
ground). The background from B Xc2 + X decay 
modes is negligibly small and ignored. 

The signal yield is extracted by maximizing a two- 
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FIG. 3; Scatter plot of AM versus AE for B° x^itt" 
candidates. The dashed hnes and solid box indicate the signal 
region. 



dimensional (2D) extended likelihood function, 
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where N is the total number of the candidate events, i 
is the event number, and Pk are the yield and proba- 
bility density function (PDF) of the component k, which 
corresponds to signal, peaking background and combina- 
torial background. The scatter plot of AM versus IS.E 
for XciT^'^ candidates is shown in Fig. 3. The 

number of events in the (Ai?, AM) fit region is 357. 

The signal PDF is modeled using a Crystal Ball (CB) 
lineshape function for Ai? and a CB lineshape func- 
tion for AAf whose shape parameters are determined 
from a signal MC sample. To take into account a 
small difference between data and MC, the shapes of 
the A_B and AM distributions are corrected accord- 
ing to calibration constants obtained from the — )■ 
J/iPK*+{K*+ K+n°) and P- -> XciK~ samples, 
respectively. In the reconstruction of J/il)K*^ 
events, we require the momentum of the i:^ to be greater 
than 0.75 GeV/c. This requirement results in a AE' dis- 
tribution similar to that for the signal events. The cali- 
bration constants for the mean and width of AE' (AA/) 
are found to be -6.23 ± 0.97 MeV (-1.16 ± 0.46 MeV) 
and 1.37 ±0.07 (1.12 ±0.05), respectively. 

The peaking background shape is modeled using a CB 
lineshape function in AM and an exponential function in 
AE. The shape parameters of the CB lineshape function 
are fixed from MC and those for the exponential function 
are floated. The shapes of combinatorial backgrounds are 
modeled by a first-order polynomial function for AE and 
a second-order polynomial function for AAf. The shape 
parameters and the number of combinatorial background 
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FIG. 4; (a) The projections in /S.E for events satisfying 
0.380 GeV/c^ < AM < 0.435 GeV/c^ and (b) the projections 
in AM for events satisfying -0.09 GeV < AS < 0.05 GeV. 
The solid curve represents the overall fit, dashed curve repre- 
sents the sum of peaking and combinatorial background, and 
dotted curve only combinatorial background component. 



events are allowed to float in the fit. 

The AE and AM distributions along with the projec- 
tions of the fit are shown in Fig. 4. The fit yields a signal 
of 40±9 B^ Xci7r° candidates. The number of peaking 
background events is 14 ± 7, which is in good agreement 
with MC expectations of 14. 

The sign ificance of the B ^ Xcitt^ signal is 4.7cr, de- 
fined as y/~2 ln(£o/'Cmax) and £max (-Co) denotes the 
maximum likelihood value (the value obtained from the 
fit when signal yield fixed to zero) . We include the effect 
of systematic uncertainties by subtracting the quadratic 
sum of the variations of the significance in smaller direc- 
tion when each fixed parameter in the fit is changed by 
±la. The branching fraction (B) for the B^ Xcitt^ 
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TABLE L Summary of systematic errors on branching frac- 
tion. 



Source 



Uncertainty (%) 



Yield uncertainty 
Tracking 

Lepton identification 
7 and tt" detection 
tt" veto 
MC statistics 

Nbb 

Secondary branching fractions 



4.5 
2.0 
3.9 
6.1 
1.6 
2.3 
1.4 
5.5 



Total 



10.8 



decay mode is calculated as follows: 
B 



e X Ngg X B 
SI 



(2) 



where TVgig is the observed signal yield, e {Ngj^) is the re- 
construction efficiency (number of B mesons in the data 
sample), and Bscc is the product of B{xci J/ipl): 
B{J/ip — > U) and B{tt^ 77). The reconstruction effi- 
ciency is determined from a signal MC sample. After a 
small correction for the muon identification requirement 
the efficiency is found to be 13.0%. We use the daughter 
branching fractions published in Ref. Equal pro- 

duction of neutral and charged B meson pairs in T(4S) 
decay is assumed. The resulting branching fraction is 



B{B° Xci7i-°) = (1.12 ± 0.25 ± 0.12) x 10"^ 



(3) 



where the first error is statistical and the second is sys- 
tematic. 



The systematic uncertainties are summarized in Ta- 
ble [H The systematic uncertainty on the signal yield 
is calculated by varying each shape parameter fixed in 
the fit by ztla, and then taking the quadratic sum of 
the deviations from the nominal value. We have checked 
for possible bias in the fitting using a MC sample; no 
significant bias was found. The systematic uncertainty 
assigned to the signal yield is 4.5%. The uncertainty on 
the tracking efficiency is estimated to be 1.0% per track, 
while that due to lepton identification is 3.9%. We also 
assign an uncertainty of 4.1% for tt" —^ 77 reconstruc- 
tion, and an uncertainty of 2.0% for the 7 detection ef- 
ficiency; these are correlated and added hnearly (6.1%). 
The systematic uncertainty due to the Xci iJ/i' and 
J/V' i~^£~ branching fractions is 5.5%. The total sys- 
tematic error is the sum of all the above uncertainties in 
quadrature. 

In summary, we report the first evidence of B^ 
Xci7r° with (657 ± 9) x 10^ BB events. The observed 
signal yield is 40 ±9 with a significance of 4.7ct including 
systematic uncertainty. The measured branching frac- 
tion is B{B° Xci7r°) = (1.12 ± 0.25 ± 0.12) x 10-^ 
which is consistent with the factorization model. 
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